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Fe IX Calculations for the Solar Dynamics Observatory 
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ABSTRACT 

New calculations of the energy levels, radiative transition rates and collisional 
excitation rates of Fe ix have been carried out using the Flexible Atomic Code, 
paying close attention to experimentally identified levels and extending existing 
calculations to higher energy levels. For lower levels, R-matrix collisional exci- 
tation rates from earlier work have been used. Significant emission is predicted 
by these calculations in the 5f-3d transitions, which will impact analysis of SDO 
AIA observations using the 94A filter. 

Subject headings: Atomic Data — Sun: UV radiation 



1. Introduction 

The launch of the Solar Dynamic Observatory (SDO) allows ob servation of t he Su n 
in unprecedented detail. The Atmospheric Imaging Assembly (AIA, iLemen et al.l (1201 ll )) 
provides multiple simultaneous images of the solar disk every 12 seconds, taken through a 
variety of narrowband filters centered on individual emission lines of interest. One such filter 
in centered on 94A, and targets both the Fe xvill (93.923A) transition from 2s^ 2p^ ^Si/2 — )■ 



2s^ 2p^ ^P3/2 at high temperatures and the Fe x (94.012A) line from the 3p^ 4s^ ^D 



5/2 



-)■ 



3p P3/2 transition in cooler plasma. The emission from these two lines occurs in very 
different temperature ranges and therefore can be distinguished if the plasma temperature 
distribution is known (jBoerner et al.ll201ll ). 



Lepson et al.l (120021 ) observed spectra of Fe Vll to Fe x using an electron beam ion 
trap (EBIT) and a grazing incidence spectrometer with resolution of ~ 300 at lOOA to 
observe lines in the 60 to 140A range. They estimated that 70% of t he emission in this band 



was unaccounted for by the ex isting atomic data in the mekal (JKaastra fc Mewd Il993 



Mewe. Kaastra fc LiedahJ 19951) database. Since then, none of the r najor atomic database s 
(e.g. CHIANTI jDere et al.ll2009h . AtomDB JPoster et al.lboilh . nist jRalchenko et al.lboilh . 
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and the successor to mekal, spex (JKaastra. Mewe fc Nieuwenhuijzenlll996[ )) have updated 
their atomic data for Fe ix to include any hues in this region. Of particular interest are two 
Fe IX lines, observed at 93.59 and 94. OTA, which both fall in the bandpass of the 94A AIA 
filter. 

These lines were identified by comparing the KBIT results with calculations using the He- 
brew University Lawrence Livermore Atomic Code (hullac, iBar Shalom. Klapish fc Oreg 
(200l|)) as belonging to the 3p^5f^ — ?■ 3p^3d^ transitions of Fe ix. The structure and col- 
lisional calculations of the current version of chianti, version 6.0.1 (and, by extension, 
AtomDB which uses the chianti data for this ion) only include configuratio ns up to the n= 4 



princ ipal quantum shell, and therefore omit these transitions. The NIST (JRalchenko et al. 



201 ll ) database does include some n=5 transitions but not these lines or levels. 



In th is work, we have carried out calculations using the Flexible Atomic Code (FAC, 



Gul (120031 )) to extend the energy level and collisional calculations to include higher energy 
levels up to the n=6 shell, merged the resulting data with the best available co llisional and 



radia tive data for lower levels, and used the coUisional-radiative code apec ( ISmith et al. 



200ll ) to model the resulting emission and therefore the effect on the AIA 94A filter flux. 



2. Method 



The NIST atomic spectra database ( JRalchenko et al.l 120111 ) lists 35 observed energy 
levels for Fe ix, including the lowest 17 energy levels and 18 others up to the 3p^ 5s^ level. 
Solar observations using the HINODE/EIS instrum ent have led to the identification of four 
more energy levels: 3p'^3d^ ^G[45^3] and 3p^4p^ ^Sq (lYoungll2009l ). 



Theoretical calculations of the Fe ix structure have been p erformed by several groups. 
The co nfigurations included in some of t hese are listed in Table [H IStorev. Zeippen &: Le Dourneul 
( 120021 ) performed superstructure (JEissner. Jones Sz Nussbaumerl Il974j ) calculations of 



the fir st 140 energy levels combined with an R-Matrix collision calculation. lAggarwal et al. 



(12006 ') revisi ted this using the General Purpose Relativistic Atomic Structure Package (GRASP, 
Dyalll (119891 )) and FAC to calculate energy levels and transition rates. They used various 
configuration combinations in order to match the first 17 energy levels as closely as possible, 
while also paying particular attention to the first 360 energy levels. They found that the 
effect of the CI between many of the configurations on the energy levels and subsequent 
oscillator strength calculations is of impo rtance, in particula r the 3s^3p^3d . This was omit- 



ted during another set of calculations by IVerma et al.l ( 12006 



code, and gave very different results compared to those of 



using the Civ3 (lHibbertlll975l) 



Aggarwal et al.l ( 12006! ) and ob- 
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Table 1. The configurations included in previous caclulations of energy levels, A- values 

and collision strengths of Fe ix. S2002 denotes the IStorey. Zeippen fc Le DourneufI (120021 ) 

data while A2006 denotes lAggarwal et al.l (120061 ) . Levels marked "CI" are included only for 

configuration interaction purposes. 



Config. 



S2002 A2006 Current Config. 



S2002 A2006 Current 



3s2 3p6 






Y 


Y 


Y 


3s2 3p3 3d3 




CI 


Y 


Y 


3s2 3p5 


3di 




Y 


Y 


Y 


3s2 3p2 3d4 








CI 


3s2 3p5 


4si 




Y 


Y 


Y 


3si 3p6 3di 




Y 


Y 


Y 


3s2 3p5 


4pi 




Y 


Y 


Y 


3si 3p6 411 






Y 


Y 


3s2 3p5 


4di 




CI 


Y 


Y 


3si 3p6 511 








CI 


3s2 3p5 


4fi 






Y 


Y 


3si 3p6 611 








CI 


3s2 3p5 


51^ 








Y 


3si 3p5 3d2 




CI 


Y 


Y 


3s2 3p5 


61^ 








Y 


3si 3p5 3di 


411 




Y 


CI 


3s2 3p4 


3d2 




Y 


Y 


Y 


3si 3p5 3di 


511 






CI 


3s2 3p4 


3di 


411 




Y 


Y 


3si 3p5 3di 


611 






CI 


3s2 3p4 


3di 


511 






CI 


3si 3p4 3d3 




CI 


Y 


CI 


3s2 3p4 


3di 


611 






CI 


3p6 3d2 




CI 


Y 


CI 


3s2 3p" 


4p^ 








CI 


3p^ 3d3 




CI 




CI 
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served energy levels. The configurations listed in Table [T] are those from their best fitting 
GRASP run. 



The energies of IStorey. Zeippen &: Le DourneufI (120021 ) agree with the observed energies 
to within 6% (the energies of the first 17 levels wer e set to match observe d values after the 
structure calculation was complete), while those of lAggarwal et al.l (120061 ) ag ree within th e 
observed energies 1 to 3%. The exception to this are the four levels from (JYoungl |2009|) , 
numbers 94, 95, 96 an d 144 in our energy ordering. These were not identified until after the 
Aggarwal et al.l (120061 ) data was created, and there is a substantial difference in both the 
energies and the level ordering caused by their introduction (see Figured]). 

In our FAC calculation, we have experimented with different configuration sets to obtain 
a good match with the observed energies in both absolute energy and energy level ordering, 
and to include the higher n shells from which the emission takes place. We have included the 
3s^3p^51 and 3s^3p^61 configurations to include the higher n emission which we are trying 
to characterize. 

We have also included many other configurations for their CI effects. The coupling 
between configurations of the same parity can affect the energy ordering of many of the levels 
during the structure calculation. In this case, most of these levels do not contribute observed 
emission lines, and their energies are much higher than those of immediate spectroscopic 
interest for which there are observed lines (> 200eV), leading them to straddle the ionization 
energy of the ion (i.e. many of them are auto- ionizing). To keep our calculation to a 
reasonable size, we have included many configurations only during the structure calculation 
for CI purposes and omitted them from the coUisional and radiative calculations: again, 
these are listed in Table [H 

Of these additional configurations, it is the inclusion of the 3s^3p^3d^ and 3s^3p^4p^ 
config urations in th e structure calculation which adds the CI necessary to bring the energies 
of the lYoung (120091 ) levels back to within 2.5% of the observed values, and also creates an 
energy ordering which matches observations for all identifie d levels. Using our mu ch larger 
FAC calculation, we obtain values comparable to those of lAggarwal et al.l (|2006[ ) for the 
lower energy levels and again within 1 to 2.6% of observed values at higher levels. This does 
le ad to significant c hange s in the energy ordering of levels when compared to the calculations 
of lAggarwal et al.l (120061 ). Our energy levels, combined with our best attempts to identify 
matching levels have been listed in Table |2l 

We note (see Figured]) that our calculated energies are without exception higher than the 
observed values. This offset scales simply with energy, resulting in calculated levels which are 
1-2.5% larger than observed values. Since the goal of these calculations is to produce useful 
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Fig. 1. — The ratio of the observed and cal culated energies for Fe ix levels from a variety 
of metho ds. Squares: supers tructure (IStorev. Zeippen fc Le Dourneuj |2002[ ): circles: 
GRASP (jAggarwal et al.l 120061): stars: this work. The levels with observed energy levels 
from HINODE observations (lYoungll2009l ) are highlighted, as is the 1% to 3% difference band 
in which all the current results fall. The dashed line indicates the correction factor used 
for all higher energy levels. A description of these levels is in Table |2l with their indexes in 
column 5. Only those levels with observed counterparts are included. 
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Table 2. The list of energy levels resulting from this work. EpAc refers to the original 
results from FAC calculations, Ecorr are the energies after correction as described in the 
text. For each level with an observed energy value, the index of this level in Figure 1 is 
given (Indfig- i) : for these levels Ecorr = Eotserved- For comparison, the grasp results of 
Aggarwal et al.l (120061 ) are listed (Eg), along with the energy order from that work (Indc) 
A star denotes a level for which a different conf i gurati on is found between our work and 





A^ 


earwal et al. (2006). 




















Index 


jj Symbol 




■t^FAC 

(eV) 


-*— 'corr 

(eV) 


Indpig.i 


Eg 
(eV) 


Indo 


1 


3p2 3pt(0,0)o 

2 2 




0.000 


0.000 





0.000 


1 


2 


3p|3p|3d|(0,|,|)^ 




50.949 


50.309 


1 


50.925 


2 


3 


3pi3p|3dl(0,|,|)^ 




51.264 


50.625 


2 


51.242 


3 


4 


3p|3p|3d| (0,1,1)2 




51.923 


51.288 


3 


51.907 


4 


94 


3pi3p|4si(0,f,i)^ 




120.922 


117.850 


17 


121.033 


94 


95 


3p2 3pl3dl3dl (0,2,1, 

"2222 


2)a 


121.280 


118.637 


18 


123.412 


97 


96 


3pl3p33 3dl3di (i,|,| 


2^5 


121.346 


118.703 


19 


123.470 


98 


97 


3pi3p|3dl3di(0,2,f, 


2)?. 


121.461 


118.761 


20 


123.521 


99 


98 


3pl3p|4sl(i,0,l)^ 




122.173 


119.247 




122.330 


95 


99 


3pi3p|3dl3di(i,f,| 


2J3 


122.577 


119.521 




124.507 


100 


100 


3p\3pi4s\(i,0,i)^ 

2 2 2 ^ ' ^ 




122.863 


119.715 


21 


122.951 


96 



Note. — Table 2 is published in its entirety in the electronic edition of ApJL. A 
portion is shown here for guidance regarding its form and content. 
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spectra for astrophysical analysis, we have adjusted the level energies after the structure 
calculation to match the observed values. For intermediate levels, the multiplier has been 
interpolated in energy and applied. The 3p^5f^ levels of primary interest in this work are of a 
higher energy than any other experimentally identified levels of Fe ix. Noting the relatively 
uniform overestimate of the energy in the structure calculations, we have scaled the energies 
for all levels above the highest observed energy level by the mean of the adjustment used 
for experimentally identif ed level s , 0.98 17. The observed energy levels are the 17 levels from 



NIST, the 4 levels from lYoungl ( l2009l ). and four levels calculated from the iLepson et al. 



(I2OO2I) measurements: (Spf 3p| 4f|)2 (E=163.2eV), (3pl 3p| 4f|)2 (E=164.9eV), (3p| 3p' 
51^7)2 (E=188.9eV), (3p\ 3p| 5^)2' (E=189.6eV). 



3 

3 
2 2 



Radiative rates were calculated by lAggarwal et al.l (120061 ) based on their structure cal- 
culation. We have initially attempted to use these radiative rates, however the significant 
disagreements in the en ergy level ordering hav e made this problematic for most levels. We 
have therefore used the lAggarwal et al.l (120061 ) values for transitions among the lowest 17 
energy levels, where the ordering is definite, and the remainder have been calculated using 
the relativistic method within FAC. 



Storey. Zeippen fc Le Dourneuj (120021 ) performed R-Matrix calculations, producing col- 
lisional excita tion rates among t he lowest 140 energy levels. Further work was performed 



using FAC by iLiang et al.l ( 120091 ). although again these did not incorporate n=5 configura- 
tions. 

A full R-matrix calculation including all of the excited levels on the n=5 level would 
be a prohibitively large calculation and is beyond the scope of this work. We have therefore 
use the distorted wave FAC collision code to calculate collision strengths between all of the 
levels of Fe ix included in our structure calculations. The advantage of this approach is 
the fast production of results and the inclusion of large numbers of levels. The downside 
is the omission of low-energy reso nance effects, which can be signi f icant. We therefore use 
the R-matrix collision strengths of IStorey. Zeippen fc Le Dourneuj (120021 ) where they exist. 
We note that the problem of matching levels between the calculations persists. For most 
levels, we have matched the levels with the same electron configuration and total angular 
momentum, J, and then paired these results in energy order. For the 3p^3d^ and 3p^4p^ 
levels there are different numbers of levels (for the forme r configuration, our calculation 
has 2 more levels for each of J=2, 3 and 4; for the latter IStorey. Zeippen fc Le Dourneui 



(120021 ) have 2 more). We have therefore combined these two configurations, which overlap 
completely in energy at the higher end of the 3p^3d^ energies (~ 130eV). 



3. Results 

We have taken the results of our FAC calculation, merged where appropriate with 
other calculations as outlined above, and used the apec code to mode l the emission from 
an op tically thin plasma with solar photospheric elemental abundances (jAnders &: Grevesse 



19891 ) in collisional ionization equilibrium. We show this result in Figure |21 for T=10^K 



plasma with the emission lines broadened by Gaussians with width 0.0 5A. Overplotted o n 



this figure is the same result from the currently version (2.0) of AtomDB (JFoster et al.ll201ll ) 



which uses the old Fe ix data. Overlaid is the effective area curve of the AIA instrument 



with the 94A filter in place, taken from iBoerner et al.l (120111 ). which can also be obtained 



by the routine aia_get_response.pro in the SolarSoftware package. It can be seen that there 
are clearly two strong lines and several weaker lines, which are also from 5f^3d transitions, 
within the region of interest. 

The weaker lines are a significant fraction of the two main 5f-3d lines, with emissivities 
of around 10-20% of the main lines. Their exact wavelengths are, however, unknown: the 
5f-3d lines were already weak in the KBIT measurements, the weaker neighboring lines were 
not distinguished from the background. The correction applied to the 5f-3d lines was a 
further ±5% compared to the simple multiplication by the 0.9817 scaling factor, in opposite 
directions. This implies that a correction at a similar level may be required in the case of 
these weaker lines. We have not further adjusted the energy levels of these weaker lines 
after the general initial scaling. Further experimental measurements of these lines would be 
valuable in estimating their effects more clearly. 

In Figure [2] we have also convolved the line emissivities for all ions with the effective 
area of the filter at each wavelength, and show the total as a function of temperature. We 
have done this for the AtomDB 2.0.1 model, which omits the 94A lines, and for an identical 
dataset but using the new Fe ix data from this work. The low temperature peak, previously 
due to Fe x, is significantly increased due to the strong Fe ix emission as part of this work, 
by over a factor of 2 at \ogiQ(Te/ K) = 5.8. 

In addition, further lines are positively identified in the Fe ix spectrum. The lines at 
134.08 and 136. 70A, identified as belonging to the 3p^4f -^ 3p^3d transition, are clearly 
identifiable in the new calculations of the spectrum. Figure [3] shows the spectrum calculated 
using the old and new Fe ix data, combined with a quiet Sun spectrum taken from the 
SDO Extreme Ultraviolet Variability Experiment MEGS. The previously unidentified lines 
at 134. 08A and 136. 70A are clearly observable. These wavelengths fall sufficiently far away 
from the 13lA filter transmission band that there is no discernable change in the estimated 
fiux in this filter due to the inclusion of these lines. 
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Fig. 2. — Left: the emissivity of Fe ix at T = lO'^K. Dotted line: from AtomDB v2.0.1; solid 
line: from AtomDB v2.0.1 with Fe IX data from this work; dashed line: the effective area of 
the 94A channel of the SDO/AIA. The AtomDB v2.0.1 data is effectively zero, and hence 
too small to be notable on the graph. Right: The total emissivitiy of elemental emission 
lines convolved with the effective area of the 94A filter. Dashed line: data from AtomDB 
v2.0; solid line: data from AtomDB v2.0 including the Fe IX data from this work. 
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Discussion 



The lack of i nformation on the Fe i x emission in the 9 4A region has been well known 
since the work of iLepson et al.l (120021 ) . iTesta et al I (l201l[ ) compared the existing atomic 
data in chianti 6.0.1 with observations of Procyon using the Chandra LETG spectrometer, 
which covers wavelengths up to ~ 170 A. They noted in these observations a significant 
missing flux in the 94A region. Since the launch of the SDO, the lack of a coUisional 
excitation calculation has led to several authors using creative means to account for this 
difference. When looking at coronal loop emission ISchmelz et al.l (J201ll ) used a scaling factor 
obtained by comparing the rel ative intensiti e s of th e known transition lines and the 5f -^ 3d 
in the spectrum observed by iLepson et al.l (120021 ). They eventually discounted this band 
from their analysis due t o a lo w count rate, so the effectiveness of this method is unknown. 
Aschwanden fc Boernerl (l201ll ). when looking at the temperature structure of coronal loops, 
derived a "correction factor", ^94, for the low temperature {\ogiQ{T/K) < 6.3) response 
function of the 94A filter, i?94, by fitting their 100 results in the other filters with ^94 a 
free parameter. They obtained ^94 = 6.7 ± 1.7 at temperatures around \ogiQ(T / K) = 6.0. 
The results from this work imply a much smaller correction factor, with ^94 ?« 2 being the 
maximum value at log;^o(^/-^) — 5-8, and ~ 1.25 at log^o(^/-^) — 6-0. There is significant 
di screpancy here: it would b e an i nteresting exercise to investigate whether the DEM analysis 
of lAschwanden fc Boernen (120111 ) can be self consistent using the new emission estimate in 
this filter band. It is also possible that there are futher lines from different ions in the filter 
which have not yet been correctly handled. 

It is difficult to estimate uncertainties in the results from this work due to the many 
overlapping sources of error. As well as the two main lines which are the focus of this 
work, there are many smaller lines of uncertain wavlength. Given the observed and initially 
calculated 5f — )■ 3d wavelengths differ by around 0.5%, adding normally distributed random 
errors to these lines of ±0.5A gives a 10% standard deviation in the resulting Fe ix emission 
in the 94A filter. This 10% fiuctuation is then approximately a 5% effect on the total 
emission in the band once Fe x is also included. Given that the estimates of uncertainty in 
distorted wave excitation collisions are usually not better than 20%, this is not expected to 
be a dominating source of error in the calculation. 

For those levels which ap p ear i n both this work and the collisional calculation of 
Storey. Zeippen fc Le Dourneuj (120021 ). comparison of the collision strengths for excitation 
from the ground state at T = lO^K between this data and our FAG results show that they 
vary by on average around 25%. This is a very approximate lower limit on the likely errors 
on the calculation of emission in the 94A band, since excluding cascades from the n=6 shell, 
this is the only way to populated the 5f upper levels in this model, and for these levels we 
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are using the distorted wave, not R-matrix, methods. 



This data has been incorporated into AtomDB, to be released fully in AtomDB v2.1.0 
which is due for release later in 2011. In the meantime, the data can be obtained from the 
AtomDB website, [www . atomdb -org 



ARF acknowledges funding from NASA ADP grant #NNX09AC71G. PT has been 
supported by contract SP02H1701R from Lockheed-Martin to SAO. 



REFERENCES 

Aggarwal, K. M., Keenan, F. P., Kato, T., and Murakamia, I., 2006, A&A, 460, 331 

Anders, E., and Grevesse, N., 1989, Geochimica et Cosmochimica Acta, 53, 197 

Aschwanden, M. J., and Boerner, P., 2011, ApJ, 732, 81 

Bar-Shalom, A., Klapisch, M., and Oreg, J. 2001, JQSRT, 71, 169 

Boerner, P., et al., 2011, to appear in Solar Physics, 



http://www.lmsal.com/sdodocs/doc?cmd=dcur §^proj_num=SDOD0063&file_type=pdf 



Dere, K. P. et al. 2009, A&A, 498, 915 

Dyall, K. G., Grant, I. P, Johnson, C. T., Parpia, F. A., and Plummer, E. P., 1989, CPC, 
55, 424 

Eissner, W., Jones, M., and Nussbaumer, H., 1974, CPC 8, 270 

Foster, A. R., et al., in prep. 

Gu, M. F., 2003, ApJ, 582, 1241 

Hibbert, A., 1975, CPC, 9, 141 

Kaastra, J. S., and Mewe, R., 1993, Legacy, 3, 6 

Kaastra, J. S., Mewe, R., and Nieuwenhuijzen, H., 1996, UV and X-ray Spectroscopy of 
Astrophysical and Laboratory Plasmas, 411 

Lemen, J., et al., 2011, Solar Physics, 10.1007/sll207-011-9776-8. 

Lepson, J. K., et al., 2002, ApJ, 578, 648 



- 12- 

Liang, G. Y., et al., 2009, ApJ, 702, 838 

Mewe, R., Kaastra, J. S., and Liedahl, D. A., 1995, Legacy, 6, 16 

Testa, P., et al, 2011, submitted to ApJ. 

Ralchenko, Yu., Krainida, A.E., Reader, J., and NIST ASD Team (2011). NIST Atomic 



Spectra Database (ver. 4.1.0), [Online]. Available: http://physics.nist.gov/asd3 [2011, 



June 15]. National Institute of Standards and Technology, Gaithersburg, MD. 
Schmelz, J. T., et al., 2011, ApJ, 731, 49 

Smith, R. K., Brickhouse, N. S., Liedahl, D. A., and Raymond, J. C., 2001, ApJL, 556, 91 
Storey, R J., Zeippen, C. J., and Le Dourneuf, M., 2002, A&A, 394, 753 
Verma, N., Jha, A. K. S., and Mohan, M., 2006, ApJS, 164, 297 
Young, R R., 2009, ApJL, 691, 77 



This preprint was prepared with the AAS E^TJtjX macros v5.2. 



-13- 



B S 3.0 

fa > 



g 2.0 



^ 4 



t5 =« 

CO G 

CO u 

s -^ 



SDO EVE quiet sun spectrum 



Fe IX lines 




AtomDB calculated spectrum 
atTe = 8xlO°K 




120 



125 130 

Wavelength (A) 



135 



140 



Fig. 3. — Top: A quiet Sun spectrum taken from the SDO EVE MEGS instrument. Bottom: 
The spectrum calculated using AtomDB 2.0.1, with (dashed) the old Fe ix data and (solid) 
the new data from this work, with all lines broadened by Gaussians with a = 0.2A. 



